Oscillations of Dirac and Majorana neutrinos in matter and a magnetic field 
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We study the evolution of massive mixed Dirac and Majorana neutrinos in matter under the 
influence of a transversal magnetic field. The analysis is based on relativistic quantum mechanics. 
We solve exactly the evolution equation for relativistic neutrinos, find the neutrino wave functions, 
and calculate the transition probability for spin-flavor oscillations. We analyze the dependence 
of the transition probability on the external flelds and compare the cases of Dirac and Majorana 
neutrinos. The evolution of Majorana particles in vacuum is also studied and correction terms to 
the standard oscillation formula are derived and discussed. As a possible application of our results 
we discuss the spin-flavor transitions in supernovae. 

PACS numbers: 14.60.Pq, 14.60.St, 03.65.Pm 

Keywords: neutrino oscillations, background matter, magnetic field, supernova neutrinos 



I. INTRODUCTION 



If massive neutrinos possess nonzero transition mag- 



that change both 



netic moments, transitions 
spin and flavor of neutrinos can happen in electromag- 
netic fields. Such transitions may be realized, e.g., in 
astrophysical environments where strong magnetic fields 
are present. Oscillations of Dirac neutrinos in an ex- 
ternal magnetic field were studied in Ref. In this 
scenario one has a transition to the sterile neutrino state 
While considering oscillations of Majorana neutrinos 
with transition magnetic moments in an external mag- 
netic field (see, e.g., Refs. one has transitions 
between active neutrino states since = {v]^Y fo^' ^ 
Majorana particle. 

Neutrino spin-flavor oscillations in solar magnetic fields 
have been earlier studied in connection to the solar neu- 
trino problem (see, e.g., Ref. [5]). It was thought that this 
neutrino oscillations channel could explain at least partly 
the deficit of the electron neutrinos in the measured solar 
neutrino flux. It is now clear that such spin-flavor oscilla- 
tions cannot play any significant role and the deficit can 
be satisfactorily explained in terms of active-active con- 
versions and the Mikheyev-Smirnov-Wolfenstein matter 
effect [Ij. 

The influence of strong magnetic fields of neutron stars 
on neutrino oscillations was studied in Refs. [7,j&.9,]- It 
was shown that spin-flavor transition may have impor- 
tant effects in the neutron star environment. In these 
investigations neutrino spin-flavor oscillations were de- 
scribed for realistic profiles of matter densities and mag- 
netic fields and the appearance of resonances in neutrino 



oscillations was examined. For further details on the neu- 
trino oscillations in electromagnetic fields as well as neu- 
trino electromagnetic properties the reader is referred to 
the recent review [lOj . 

In this paper we shall return to the question of the spin- 
flavor oscillations using an approach that differs from 
the one usually followed. It is usual to describe neu- 
trino oscillations, including the spin-flavor transitions, on 
the basis of the quantum mechanical evolution equation. 
Instead of this Schrodinger picture, we will apply rela- 
tivistic quantum mechanical picture based on the Dirac 
theory. In this approach, we will study spin-flavor oscil- 
lations of Dirac and Majorana neutrinos in matter and 
in an external magnetic field, extending our earlier use 
of the method [ll|,1il [H, [13 to this new problem. We 
should also mention that the majority of the previous 
studies of neutrino spin-flavor oscillations are restricted 
to the case of Majorana neutrinos. We will investigate 
both the Dirac and Majorana cases as the nature of neu- 
trinos is still an open question 
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The plan of this paper is as follows. We will start by 
writing down the relativistic wave equations that take 
into account the background matter and external mag- 
netic field. We formulate the initial condition problem for 
these systems (see also Refs. [H [H, [H We then 
derive for both Dirac and Majorana neutrinos a Hamil- 
tonian analogous to that of the standard quantum me- 
chanical approach and solve exactly the resulting evolu- 
tion equation for relativistic neutrinos. We then analyze 
the behavior of the transition probability for Dirac and 
Majorana neutrinos at various magnetic field strengths 
and matter densities. Our results will be summarized 
and their applications to some astrophysical situations 
discussed in Sec. IVl 
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II. EVOLUTION OF DIRAC NEUTRINOS IN 
MATTER AND TRANSVERSAL MAGNETIC 
FIELD 

Let us study the evolution of two neutrino flavor states 
h'x, X ~ a,f3, in a nonnioving and unpolarized matter 
under the influence of an external magnetic field. We as- 
sume that the mass eigenstates ipa, a, — 1,2, are related 
to the flavor eigenstate neutrinos through the transfor- 
mation 



a=l,2 



cos 9 — sin 9 



(2.1) 



where is a mixing angle. In this section we shall assume 
that the mass eigenstates are Dirac particles. We set the 
following initial conditions (see also Refs. [11, 12. 13, 14]): 



i^„(r,0)=0, iy0ir,O) 



(0) „ikr 



(2.2) 



(0)T 



where k = (k, 0, 0) is the initial momentum and — 
(1/2)(1, — 1, — 1, 1) (see below). This corresponds to a 
relativistic neutrino of the flavor /?, with its spin directed 
oppositely to the particle momentum, i.e. a left-handed 
neutrino v^. The system is taken not to contain the other 
neutrino flavor Va initially. 

Note that from the physical point of view it would be 
more realistic to choose a localized in space initial con- 
dition rather than that in Eq. (|2.2[) . It is well known, 
however, that from the point of view of oscillations the 
plane wave and wave pac ket approaches in practice lead 
to equivalent results [16| . Moreover the initial condi- 
tion problem (|2.2[) for the cases of flavor and spin-flavor 
oscillations in various external fields was solved in our 
previous publications [H!, [l2, 13, 14]. 

The Dirac equation for the mass eig enstates wave func- 
tions is of the form (see Refs. [isl fl^). 

i^a = Hai^a + {Vm + l^s)^fc, a ^ b, (2.3) 

where Tia — (cip) + -I- (70(1 — 75) /2 is the mass- 
diagonal part of the Hamiltonian, ma are the masses asso- 
ciated with the states V'a, and a. = 7^7, /3 = 7", and S — 
7^7°7 are Dirac matrices. The term Vm = g(l — 75)/2 
describes the interaction of neutrinos with particles in 
the background matter. The matrix (gab), given below, 
is not diagonal in the mass eigenstate basis. The term 
Vb = —/ii? /3S3 is the energy operator associated with the 
interactions of neutrinos with the magnetic field, and the 
magnetic moment matrix (fiab) is, like {gab), in general 
nondiagonal. The magnetic field is taken to be transver- 
sal with respect to the initial neutrino momentum, i.e., 
B = (0,0, i?). The Vm and Vb terms are responsible 
for the possible mixing between different neutrino mass 
eigenstates. 

The matrix {gab) that describes the interactions of neu- 
trinos with matter is in the mass eigenstates basis of the 



form 

gab = uljxUxb 

„2 a I J-.. „;„2 



(2.4) 



fa cos^ 9 + sin 9 — sin 9 cos 9 [fa — ffj] 
- sin 9cos9[fa~ fp] fa sin^ 9 + fp cos^ 9 

ja and .912 = 521 = 



We will denote gaa = 
identify the flavor a as 
= Ue, the effective potentials f\, A 
(see Ref. [13): 



J. If we 
or and the flavor f3 as 
a, /3, are given by 



f=e,p,n 

qf^H4L-2Q^fhm'9w) 



if -2Q^^^ sin' 9w + Sef), 



rif) 



if), 



(2.5) 



( f) 

where Uf is the number density, /g^ the third isospin 
component and Q'^-^'^ the electric charge of the back- 
ground particle of the type /, 9w is the weak mixing 
angle, and Gp is the Fermi constant. It is assumed in 
Eq. (|2.5p that matter consists of electrons, protons, and 
neutrons and that it is unpolarized and at rest. 

As to the interactions of neutrinos with the magnetic 
field, we have assumed when writing Eq. (|2.3p that the 
magnetic moment matrix in the mass eigenstates basis is 
antidiagonal, i.e., fj,aa — and /ii2 — /i2i = 7^ 0. Mag- 
netic moment matrices of this type were studied in our 
previous works [H, [13]. The situation with the values 
of Dirac neutrino magnetic moments is disputable p^ : 
however, the existence of large off-diagonal magnetic mo- 
ments of Dirac neutrinos is not excluded. Let us note 
that the solution to the Dirac-Pauli equation for a neu- 
trino with Haa 7^ propagating in an arbitrary moving 
and polarized medium was recently obtained in Ref. 0| . 
However in that work the case of only one neutrino flavor 
was studied. 

The general solution of Eq. (|2.3p has the form [3, ]T^, 



V'a(r, = exp (-ic/at/2) 



d3p 



(27r)3/2 

X J2 [ai^\t)u(Pe^pi-iEiOt) 
C=±i 

+ bi<\t)vi^^exp{+iEi<h) 
where the energies are given by (see, e.g., Ref. [2(|) 



(2.6) 



In the relativistic limit one has 



(2.7) 



(2.8) 



To obtain Eq. (12. 8p we neglect the term ^ g^ compared to 
the neutrino mass squared . For the situation of neu- 
trino propagation in the expanding envelope formed after 
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the supernova explosion, which is studied below, the ef- 
fective potentials are ga ^ 10~^^ — 10"^^ eV. Accounting 
for the neutrino mass in the eV range, the approximation 
made is always valid. 

The expressions for the basis spinors u"P and wi^' can 
be found in Ref. [lOj- When one studies the situation 
where the initial momentum is very large, k ^ m^, one 
can neglect the neutrino mass dependence of the spinors 
and present them in the form 




(2.9) 



Note that the coefficients ai^\t) and bi^\t) in Eq. (HH) 
are functions of time due to the presence of the non- 
diagonal interactions Vm and Vb- Our main goal is to 
determine these coefficients so that both Eq. (|2.3p and 
the initial condition (12.21) are satisfied. 



Using Eqs. (|2.3p and (|2.6p and the orthonormality of 
the basis spinors (|2.9p we arrive to the following ordinary 
differential equations for the coefficients ai^^ and [H, 

iai^) = exp [i{ga ~ gb)t/2] exp {+iE[^h)u^^^\V„, + Vb) 



C'=±i 



i&l« 



: exp 



[i{ga - gb)t/2] exp i-iE^pt)v^'^'>\Vm + Vb) 



X J2 [ai'''^u^'^")exp{-iEl^'^'h) 
C'=±i 

+ b[^'K^'^'^exp{+iEl'^'h) . 



(2.10) 



One easily sees that (u \Vm\u ) ~ g, {v^\Vm\v^) — g, 
{u^\Vb\u^) = -M^ and {v ^\Vb\ u^) = -fiB. All the 
other scalar products in Eq. (I2.10p vanish. 

Let us introduce a four-component wave function 
defined as ^' = (oj^, Oj^, a|j^, a^). Then we can rewrite 
Eq. (|2.10p in the form of a Schrodinger equation: 



d*' 

i— = H'-if', H' = 
dt 






ge"^* 





-fj.Be''^'*\ 


ge-''^'* 













-fiBe''^^* 
















/ 



(2.11) 



where [see also Eq. p. 8 



uji = - E. 



+ , 91 - 92 



oj' = E^ - £'2 



2$ -I- 51, u;2^E+-E- 
91 - 92 



91 - 92 



2$ - 



2$ + 2(.gi-(72), 



(2.12) 



and $ = 5m? /{Ak) = {ml — m\)/{Ak) is the vacuum oscillation phase. 
Let us next do the following transformation to the wave function: 



U =diag |e 



i(*+3gi/4-g2/4)t „-i(*-3g2 /4+si /4)t „i(*-Si/4-32/4)t „-i(*+3i/4+32/4) 



, e 



This takes the Schrodinger Eq. (j2.1ip into the form 

(Ml . i • 

i—=H^, H=WH'U-iU^U = 
dt 



/$ + 3gi/4 - 52/4 g -^lB 

9 -$ + 3.92/4-31/4 -iiB 

-fiB $-(gi+g2)/4 

V -/iS -$-(gi+.92)/4 



(2.13) 



r 



It is in order at this stage to compare the time evo- 



lution Eq. (j2.13p . which we obtained to in our approach 
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based on Dirac equation, with the one obtained in the 
conventional quantum mechanical formalism. To this end 
let us introduce a quantum mechanical wave function 
^QM — (V"! ! V'2 ' V'l^j V'f') the neutrino mass eigen- 
states. Here V'a ^'^^ one-component objects. This wave 
function obeys the Schrodinger equation with the follow- 
ing effective Hamiltonian (see, e.g., Ref. 0]): 



QM 



^/Ci+.gi g -[iB\ 

9 K.2+92 -f^B 
-nB Ki 
-HB /C2 



(2.14) 



where /Co = -y/P 



k + m1_/{2k) is the kinetic 



energy of a massive neutrino. It is easy to see that this 
matrix leads to the same dynamics as the Hamiltonian in 
Eq. (|2.13p . Indeed, without changing the dynamics one 
can make the following replacement Hqm — > Hqm — 
tr(iJQjv/)/4 • /, where / is the 4x4 unit matrix. The 
resulting Hamiltonian is exactly the same Hamiltonian 
we have derived in our approach. 

To describe the evolution of the system (I2.13|) in a 
general case, one has to solve a secular equation which is 
the fourth-order algebraic equation (see also below). Al- 
though one can express the solution to such an equation 
in radicals, its actual form appears to be rather cum- 
bersome for arbitrary parameters. If we, however, con- 
sider the case of a neutrino propagating in the electri- 
cally neutral isoscalar matter, i.e. Ue = Up and Up = ri„, 
a reasonable solution is possible to find. We will demon- 
strate later that it corresponds to a realistic physical 
situation. As one can infer from Eq. (|2.5p for the case 
of the — > oscillations channel, in a medium with 
this property one has the effective potentials fa = fp, — 
= — Gprt/v^ and fp = fe = Ve = Gpn/V^, where 
n = He = Hp = Un- Using Eq. (12. 4|) we obtain that 
gi — —92 — 9o, where go = —T^ cos 20, g = y sin 20, and 
V = Gpn/v^. 

Let us point out that background matter with these 
properties may well exist in some astrophysical environ- 
ments. The matter profile of presupernovae is poorly 
known, and a variety of presupernova models with differ- 
ent profiles exist in the literature (see, e.g. Ref. [2]|). 
Nevertheless, electrically neutral isoscalar matter may 
well exist in the inner parts of presupernovae consisting 
of elements heavier than h ydr ogen. Indeed, for example, 
the model W02Z in Ref. [Ill predicts that in a 15Mq 
presupernova one has Ye = n-e/(np -\- n„) = 0.5 in the 
0-|-Ne-|-Mg layer, between the Si-|-0 and He layers, in 
the radius range (O.OO7-O.2)i?0. 

In this kind of background matter the effective Hamil- 
tonian in Eq. (|2.13p is replaced by 



H = 



r 


+ .90 


9 


-fiB 




9 


-($ + go) -f^B 










-fiB $ 





-fiB 





-$ 



(2.15) 



We now look for the stationary solutions of the 
Schrodinger equation with this Hamiltonian. After a 



straightforward calculation one finds 

C=±i 

+ (Vc<»V^)exp{i£ct)\ *o, 



(2.16) 



where we have denoted 



£± =- v/2V^2 + 4(^5)2 + 4$2 _ 4$^ cos 29 ± 2VR, 



R ^\/iV - 2$ cos 29f + AifiBY' 



(2.17) 



The vectors U± and V± are the eigenvectors correspond- 
ing to the energy eigenvalues £± and — f±, respectively. 
They are given by (C = ±) 



sm2e{£(; 



where 



V + CR 



piB sin 2Q 
\~tiBZj{£^ + $) 

- sin 2e{£i^ - $) 

fiBZji£^ + $) 
-HB sin 20 



- £q cos 29, 
{pBf 



(2.18) 



{£c + <^f 
sin2(26l) [{^lBf 



(2.19) 



It should be noted that Eq. (|2.16p is the general solution 
of Eq. (|2.13p satisfying the initial condition ^'(0) = »&o. 

Note that we received the solution (|2.16p - (|2.19p of the 
evolution Eq. (|2.13p under some assumptions on the ex- 
ternal fields such as isoscalar matter with constant den- 
sity and constant magnetic field. We mentioned above 
that our method is equivalent to the quantum mechani- 
cal description of neutrino oscillations |3|] which is valid 
for a more general case of coordinate dependent exter- 
nal fields. The advantage of our approach consists in the 
fact that one can derive neutrinos' wave functions for ar- 
bitrary initial momenta, as it was made in Refs. pUfl^. 
and study the propagation of low-energy neutrinos. The 
assumption of constant matter density and magnetic field 
is quite realistic for certain astrophysical environments 
like a shock wave propagating inside an expanding enve- 
lope after a supernova explosion (see also Sec. IV)) . 

Consistently with Eqs. (|2.ip and (|2.2p . we take the ini- 
tial wave function *(0) = vJ/q in Eq. (|2.16p as *J = 
(V^2>V'f:V'|^) (sin 0, cos 6*, 0,0). Using Eqs. ^f^- 
(I2.19P one finds the components of the quantum mechani- 
cal wave function corresponding to the right-handed neu- 
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trinos to be of the form 



sin 6* sin 26* Z. 



+ 



+ {+ 



Nl 



sm ( 



-i£+t 



sin2(26')(f+ - $)e 
Zi 



cos 61 sin 26* Z_ 



+ 



i£+t 



+ { + 



(2.20) 



where the {+ — > — } stand for the terms similar to 
the terms preceding each of them but with all quan- 
tities with a subscript + replaced with corresponding 
quantities with a subscript — . The wave function of 
the right-handed neutrino of the flavor a, i^^, can be 
written with help of Eqs. ((27T|) and ([2?20l) as iy^{t) = 
cosei}f {t) -sinO^/jfit). 

The probability for the transition vj^ is obtained 

as the square of the quantum mechanical wave function 
v^. One obtains 

P^L^^R (t) = [C+ cos{£+t) + C- cos{£-t)f 

+ [S+ sm{£+t) + S- sin(£_i)]2, (2.21) 



-sin2(26')(£(; -$) 




sin^(26')(£^ -$) 



(2.22) 



As a consistency check, one easily finds from Eq. (|2.22p 
that C+ -|- C_ = as required for assuring P(0) — 0. 



In the following we will limit our considerations to the 
case « corresponding to the situations where 
the effect of the interactions of neutrinos with matter 
{V) is small compared with that of the magnetic inter- 
actions (fJ-B) or the vacuum contribution ($) or both 
[see Eq. (|2.17|) ]. Note that in this case one can ana- 
lyze the exact oscillation probability (j2.21|) analytically, 
which would be practically impossible in more general 
situations. 

In the case sa one can present the transition 
probability in Eq. (|2.2ip in the following form: 

P{t) = Po{t) + Pc{t) cos(2m) + Ps{t) sin(2f7t), (2.23) 
where 

-Po(t) =^ [5"+ + 81+ 2S+S- cos{2Snt) 

-4C+C- sin^iSnt)], 
Pc{t) = - i [{Si + S^_) cos{26nt) + 2S+S- 



-4C+C_ sin2(,5m)], 
Ps(t) ^^{Sl~Sl)smi2Snt), 



and 



sn 



(2.24) 



(2.25) 



As one can infer from these expressions, the transition 
probability P{t) is a rapidly oscillating function, with the 
frequency Q, enveloped from up and down by the slowly 
varying functions Pu,d — Po +^ \J Pc + Ps j respectively. 

The behavior of the transition probability for various 
matter densities p and the values of /iB and for a fixed 
neutrino energy of i? = 10 MeV and squared mass differ- 
ence of 5m? = 8 X 10^^ eV^ is fllustrated in Figs. [T]|3l 

As these plots show, at low matter densities the en- 
velope functions give, at each propagation distance, the 
range of the possible values of the oscillation probability. 
At greater matter densities, where the probability oscil- 
lates less intensively, the envelope functions are not that 
useful in analyzing the physical situation. 

One can find the maximum value of the upper envelope 
function, which is also the upper bound for the transition 
probability, given as 



p(max) ^ I ('5+ - 5"-)^, \iB<B', 

\C+C.{Sl~Slf[C+C.{Sl + Sl) + {C+C-f + [S+S-f]-\ ifB>B', ^ ' 



where the value B' is the solution of the transcendent maximum values of the averaged transition probablility 
algebraic equation, C+C- = 5+5-. The corresponding 
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FIG. 1: (a)-(c) The transition probability versus the distance passed by a neutrino beam in matter with the density p = 10 g/cc; 
(a) nB = 5x 10"" eV, (b) = 1.1 x 10"^^ eV, (c) fiB = 5 x 10"^^ eV. We take that = lOMeV, W = 8 x lO"'^ eV^ and 
8 = 0.6, which is quite close to the solar neutrinos' oscillations parameters. The black line is the function P{x), the blue and 
green lines are the envelope functions P„_tj(a;), and the red line is the averaged transition probability Po{x). (d) The dependence 
of the maximal values of the functions P{x) and Po{x), blue and red lines, respectively, on the magnetic energy for the 
given density. 



Po{x) are given by 

p^(max) ^ 1 [(^^5_)2 _ 4C+C.], (2.27) 

for arbitrary values of B. The values of these maxima de- 
pend on the size of the quantity /iB. These dependencies 
are plotted in Figs. md)-[3l^d). In the case of rapid oscil- 
lations the physically relevant quantities, rather than the 
maxima, are the averaged values of the transition prob- 
ability, which are also plotted in these figures. 

As Figs. [lld)-[31d) show, the interplay of the matter 
effect and the magnetic interaction can lead, for a given 
magnetic moment /i, to an enhanced spin-flavor transi- 
tion if the magnetic field B has a suitable strength rela- 
tive to the density of matter p. In our numerical examples 
this occurs at /iSmax = 1.1 x 10^^^ cV for p = 10 g/cc, 
at pB„^ax = 6.6 X 10^^^ eV for p — 50 g/cc, and at 
A^-Smax = 8 X 10"-^^ eV for p — 100 g/cc. For these values 
of pB both the maxima and the average of the transi- 
tion probability become considerably larger than for any 



other values of pB. Figures [T]^b)-[3l^b) correspond to the 
situation of maximal enhancement, whereas Figs. [Ha)- 
^sl) and Figs.[T]Jc)-[2Ic) illustrate the situation above and 
below the optimal strength i?max of the magnetic field. 

It is noteworthy that the enhanced transition proba- 
bility is achieved towards the lower end of the pB region 
where substantial transitions all occur, that is, at rela- 
tively moderate magnetic fields. At larger values of pB 
the maximum of the transition probability approaches to- 
wards cos^(20). Indeed, if pB ^ max($, V), the transi- 
tion probability can be written in the form (see Ref. p[3j ) 
P{t) = cos2(20) sin^ {pBt). It was found in Ref. ^ that 
neutrino spin-fiavor oscillations can be enhanced in a very 
strong magnetic field, with the transition probability be- 
ing practically equal to unity. This phenomenon can be 
realized only for Dirac neutrinos with small off-diagonal 
magnetic moments and small mixing angle. As we can 
see from Figs. [lId)-|3Ud) the situation is completely dif- 
ferent for big off-diagonal magnetic moments. 

One should notice that for long propagation distances 



7 




0.2 0.4 0.6 0.8 1 



4 

X, 10 km 
(a) 




3 

X, 10 km 




0.8 
0.7 
0.6 
0.5 
S 0.4 
0.3 
0.2 
0.1 



0.2 



0.4 



0.6 



^B, 10" 
(d) 



-11 



eV 



0.8 



FIG. 2: The same as in Fig. [T] for the density p = 50g/cc; (a) = 3.5 x 10"^^ eV, (b) fiB = 6.6 x 10"" eV, (c) fiB = 
5 X lO^^^eV. 



consisting of several oscillation periods of the envelope 
functions, the enhancement effect would diminish con- 
siderably due to averaging. In the numerical examples 
presented in Figs. [T][3] the period of the envelope function 
is of the order of 10'^ — 10"' km, which is a typical size of 
a shock wave with the matter densities we have used in 
the plots (see, e.g., Ref. [23i]). Thus the enhanced spin- 
flavor transition could take place when neutrinos traverse 
a shock wave. 

Let us recall that the above analysis was made by as- 
suming neutrinos to be Dirac particles. We will see below 
(see Sec. lIVI) that the corresponding results are quite dif- 
ferent in the case of Majorana neutrinos. 



III. EVOLUTION OF MAJORANA NEUTRINOS 
IN VACUUM 

We now move to consider Majorana neutrinos, and 
we shall start by applying our formalism to the or- 
dinary vacuum oscillation of two Majorana neutrinos. 
The left-handed chirality component of a flavor neutrino 
i/^ = (1/2)(1 — 7^)i^A is related to the wave functions of 



Majorana neutrino states through 

a 

where A = a,/? is the flavor index and rja, a = 1,2, 
correspond to a Majorana particle with a deflnite mass 
nia- In the simplest case the mixing of the flavor states 
arises purely from Majorana mass terms between the left- 
handed neutrinos, and then the mixing matrix U\a is a 
2x2 and unitary matrix, i.e., a = 1, 2 and, assuming no 
CP violation, it can be parametrized in the same way as 
in Eq. (lO) . 

We study the evolution of this system with the follow- 
ing initial condition [see also Eq. (|2.2p ]: 

i.L(r,0) = 0, i/L(^.,o) = ^.(°)e*^ (3.2) 

where k — (0, 0, k) is the initial momentum and i^^ — 
(0,1). The initial state is thus a left-handed neutrino 
of flavor (3 propagating along the z-axis to the positive 
direction. 

As both the left-handed state i/^ and Majorana state 
rja have two degrees of freedom, we will describe them in 
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FIG. 3: The same as in Fig. [T]for the density p = lOOg/cc; (a) /iS = 4xlO~^^eV, (b) fiB = 8xlO~^^eV, (c) ^lB = 5xlO~^^eV. 



the following by using two-component Weyl spinors. The 
Weyl spinor of a free Majorana particle obeys the wave 
equation (see, e.g., [23|). 

i?'ya + {<Tp)lla + i"laO-2?7l = 0. (3.3) 

Thegeneral solution of this equation can be presented 
as [23 



?7a(r,i) 



^!P_pipr 

(27r)3/2 



C=±i 

+ ai«*(-p)«K)(-p)e'^"*], (3.4) 

where Ea = -y^m^ + |pp. The basis spinors u^a'' and ui^"* 
have the form 

^a(P) ^ "1°, <(p) =AaW'-, (3.5) 

-f^a + iPl 



«a (P) =^a- 



'Ea + \p 

where w± are helicity amplitudes given by (26j 

^e-'''^/2cos(z9/2)\ 
e»<^/2sin(i9/2) y ' 

-e-*"^/2 sin(z9/2)\ 
e''-^/^ cos(z9/2) j ' 



■w+, 



(3.6) 



the angles and giving the direction of the momentum 
of the particle, p = |p|(sini9cos0, siniJsin^, cosi9). The 
normalization factor Aa in Eq. p.5|) can be chosen as 



= 1 - 



(Ea + \p\y 



(3.7) 



Let us mention the following properties of the helicity 
amplitudes w±: 



{crp)w± = ±|p|w±, icr2W± = TWzp, 
w±i~p) = iw^{p), 



(3.8) 



(w+ (g) wl) — (w- (g) W^) = 1(72, 



which can be immediately obtained from Eq. (|3.6p and 
which are useful in deriving the results given below. 

The time-independent coefficients aj(p) in Eq. 



9 



have the following form [25| : 
1 



4(P) = 



(27r)3/2 
im 



Oa (P) 



Ea + \P 
1 



(27r)3/2 
im 



kHpW^Hp) 

^^^\-p)u-{-p) 



(3.9) 



where yyi''-' (p) is the Fourier transform of the initial wave 
function ija, 



r?i°)(p)= J d3pe-'P--^(0)(r). 



Using Eqs. (|3.4p - (|3.9[) we then obtain the following ex- 
pression for the wave function for the neutrino mass 
eigenstates: 



Va{r, t) 



d?p 
(27r)3 











[1 










i^a + |p|J 



X w_ (Xi w: 



\E„t 



Ea + \p 
X {w+®wl^ Wa\p) 



Ea + iPl 



(3.10) 



From Eqs. (|3.8p and (|3.10p it follows that a mass 
eigenstate particle initially in the left-polarized state 
ri^\r) ~ w_(k)e*'" is described at later times by 



{ h' - 


ma 




Ea + \k\\ 



2i- 



Ea + |k| 



sin(£;at)e-'*'''w+(k) }. (3.11) 



Let us notice that the second term in Eq. p. lip describes 
an antineutrino state. Indeed the spinor w+(k) satis- 
fies the relation, (crk)w+(k) = |k|u'+(k), see Eq. (13. 8p . 
Therefore it corresponds to an antiparticle, see Ref. [27| . 
This term is responsible for the neutrino-to-antineutrino 
flavor state transition vjj ^ {^^aT- 

According to Eq. p.ip and (|2.ip , the wave function of 
the left-handed neutrino of flavor a is = cosf^ry}" — 
sin 6*7/^. From Eqs. (PTT|) and ([XTT]) it then follows that 



the probability of the transition 
given by 



v^, m vacuum is 



■cos(|k|i) sin($i) 



4|k|2 

X [ml sm{Eit) - ml sin(i?2i)] 



O 



ma 

Ikl 



(3.12) 



The leading term reproduces the familiar oscillation for- 
mula of Pontecorvo describing the transitions between ac- 
i/^. The corrections to Pontecorvo's 
^ ^, and in our pre- 
vious papers [ill, [Hj Elj [l3l we derived the analogous 
corrections for Dirac neutrinos both in vacuum and in 
various external flelds. 

Analogously we can calculate the transition probability 
for the process vj^ {^aT using the second term in 
Eq. (I31T]) . 



five neutrinos 
formula were obtained first in Ref. [2. 



X [mi sm{Eit) — m2 sin(i?2i)]^ 



(3.13) 



Note that the next-to-leading term in Eq. (|3.12p and lead- 
ing term in Eq. (j3.13p have the same order of magnitude 
~m2/|k|2. 

Before moving to consider Majorana neutrinos in mag- 
netic fields we make a general comment concerning the 
validity of our approach based on relativistic classical 
field theory. It has been stated [2] that the dynamics of 
massive Majorana fields cannot be described within the 
classical field theory approach due to the fact that the 
mass term of the Lagrangian, r]'^ia2Ti, vanishes when rj is 
represented as a c-number function. Note that Eq. (j3.3p 
is a direct consequence of the Dirac equation if we sug- 
gest that the four-component wave function satisfies the 
Majorana condition. Therefore a solution to Eq. p.3p . 
i.e., wave functions and energy levels, in principle does 
not depend on the existence of a Lagrangian resulting 
in this equation. The wave equations describing ele- 
mentary particles should follow from the quantum field 
theory principles. However quite often these quantum 
equations allow classical solutions (see Ref. 29] for many 
interestin g examp les ] . We have also demonstrated in 
Refs. [mm, O that oscillations of Dirac neutri- 
nos in vacuum and various external fields can be de- 
scribed in the framework of the classical field theory. The 
main result of this section was to show that the quantum 
Eq. (|3.3p for massive Majorana particles can be solved 
[see Eq. (|3.1ip ] in the framework of the classical field 
theory as well. 
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IV. EVOLUTION OF MAJORANA NEUTRINOS 
IN MATTER AND TRANSVERSAL MAGNETIC 
FIELD 

For describing the evolution of two Majorana mass 
eigenstates in matter under the influence of an external 
magnetic field, the wave Eq. (|3.3p is to be modified to 
the following form: 

. . , / 9a \ , . * 9 

- ifi{(TB)a2eabVb =0: a^b, (4.1) 

where eab = i(o'2)a6, and and g were defined in con- 
nection to Eq. (j2.3p . Note that Eq. (j4.1|) can be formally 
derived from Eq. (j2.3p if one neglects vector current in- 
teractions, i.e., replace (1 — 7^)/2 with —7^/2, and takes 
into account the fact that the magnetic moment matrix of 
Majorana neutrinos is antisymmetric (see, e.g., Ref. [soj). 
We will apply the same initial condition (|3.2p as in the 
vacuum case. It should be mentioned that the evolution 
of Majorana neutrinos in matter and in a magnetic field 
has been previously discussed in Ref. [3ll |. 

The general solution of Eq. (|4.ip can be expressed in 
the following form: 



as 



?7a(r,t) 



(27r)3/2 

X K'^Hp,t)4^Hp)exp(-ii?(^)i) 

C=±i 

-fa(C)*(-p,i)z;«)(_p)exp(ii?(C)t)], (4.2) 

where the energy levels are given in Eq. (12. 7|) (see 
Ref. [i^])- The basis spinors in Eq. (|4.2|1 can be chosen 



«:(p) = -K 



w+, 



E+ + i\p\-ga/2) 

Ea + (|p| +9a/2) 

Kip) =Kw^, v+{p) = A+W_, 
where the normalization factors Ai''"' , C = i are given by 

(aF)-^ = i- 



[Ea + {\P\-Cga/2W 

Let us consider the propagation of Majorana neutri- 
nos in the transversal magnetic field. Using a similar 
technique as in the Dirac case in Sec. [Til and assuming 
k ^ TTia, we end up with the following ordinary differen- 
tial equations for the coefficients ai''''. 



d\&' 



( 

ge-'^-72 



.ge"*-72 








(4.3) 



-fiBe'"-* -ge'*+72 

\nBe^"'+* -.ge-''^+72 / 
where 'I' — (a^, aj^, a^, a^) and 



(5± ^Ef - E^ w 2$ =F 
cr± =Ef - « 2$ ± 



2 ^ 

51 +92 



By making the matrix transformation 



=U^, U = diag |ei(*+9i/2)t^ e-i(*-92/2)t^ ei(*-9i/2)t^ g-i(*+3./2)t| ^ (4 4) 

we can recast Eq. (I4.3P into the form 



(4.5) 

^lB 

V 



5/2 





fiB 


$ + 32/2 


-fiB 







$-51/2 


-g/2 





-5/2 


-$-52/2 



Let us note that the analogous effective Hamiltonian 
has been used in describing the spin-flavor oscillations 
of Majorana neutrinos within the quantum mechani- 
cal approach (see, e.g., Ref. Q) if we use the basis 

*5M = (#,V'2S[V'^']M^2^n. 

Note that the consistent derivation of the master 
Eq. (|4.ip should be done in the framework of the quan- 
tum field theory (see, e.g., Ref. |2|), supposing that the 



spinors r]a are expressed via anticommuting operators. 
This quantum field theory treatment is important to ex- 
plain the asymmetry of the magnetic moment matrix. 
However, it is possible to see that the main Eq. (|4.ip 
can also be reduced to the standard Schrodinger evolu- 
tion Eq. (|4.5p for neutrino spin-flavor oscillations if we 
suppose that the wave functions rja are c-number ob- 
jects. That is why one can again conclude that classical 
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and quantum field theory methods for studying Majorana 
neutrinos' propagation in external fields are equivalent. 

Let us again consider the situation when rig = Up ^ 
Un = n, which results in gi = —32- In this case the 
eigenvalues of the Hamiltonian (|4.5|) A = ±f ± are given 
by 



1 



R cos 20)2 + (^B)2, 



(4.6) 



where V = Gprt/VS as in SecHT] The time evolution of 
the wave function is described by the formula, 



C=±i 



(V^ (g) F/) exp {i£(;t) 



(4.7) 



where and V"^ are the eigenvectors of the Hamilto- 
nian ()4.5|) . given as 





(4.8) 



where 



< 



+ V"cos26'/2 
V sin 26 



l + ^V^^sin^(2.) 
2(^5)2(5^ + $) 



2(£:(; + $ + T/cos26i/2) 

V, 



=-[2£i^{£Q + - 1/72 + $!/ cos 26*] cos 261 

+ C,RV{£q + $ - F cos 26*72) . (4.9) 

The normalization coefhcient Nq in Eq. (|4.8p is given by 
iVc = ^l + xl+yl + zl 

Proceeding along the same lines as in Sec. |TT1 we ob- 
tain from Eqs. (PTTjl and (li?7)l - (|Tg|) the probabihty of the 
process — > as, 

P^L_^R(t) =[C+ cos{£+t) + C_ cos(£_<)]2 

[5"+ sin(£+t) -I- S- sm{£-t)f, (4.10) 

where 

Q = - -^[sin26'(a;^ + yi^z^) + cos2e{yQ - Xi^ZQ)], 



(4.11) 



C 



Consistently with Eq. (|3.2p . we have taken the initial 
wave function as 



(sin 6*, cos 61, 0,0). 



(4.12) 



With help of Eqs. (|4.9p and (|4.1ip it is easy to check that 
C+ C_ = guaranteeing P(0) = 0. 

Note that formally Eq. (|4.10p corresponds to the tran- 
sitions — > v^. However, virtually it describes os- 
cillations between active neutrinos 



since 



,L\c 



for Majorana particles. 



As in the previous case of Eq. (|2.2ip . Eq. (|4.10p can 
be treated analytically for relatively small values of the 
effective potential V. The ensuing envelope functions 
Pu,d = -Po ± y/ Pq + Ps depend on the coefhcients Cq 
and 5*,; in the same way as in Eq. (|2.24p . The transition 
probabilities at various values of the matter density and 
the magnetic field are presented in Fig. [H 

Despite the formal similarity between Dirac and Majo- 
rana transition probabilities (see Eqs. p.22p and (|4.1ip ) 
the actual dynamics is quite different in these two cases, 
as one can see by comparing Figs. [Hd)-[3fd) and Fig. |4l 
panels (b), (d) and (f). In particular, in the Majorana 
case Pu™'^^^ = 4|C+C_| for arbitrary to be compared 
with Eq. (|2.26p . while the function p^j™'^'^^ has the same 
form as in the Dirac case given in Eq. I|2.27p . In contrast 
with the Dirac case, the averaged transition probability 
does not achieve its maximal value at some moderate 
magnetic field i?max value, but both P™'^^ and p^'^^ are 
monotonically increasing functions of the strength of the 
magnetic field with 1 and 1/2 as their asymptotic values, 
respectively. We can understand this behavior when we 
recall that, at ^ max($,y), the effective Hamilto- 
nian Eq. (|4.5p becomes 



i7^ = 



/o r 
00-10 
0-100 

Vi 



(4.13) 



The Schrodinger equation with the effective Hamilto- 
nian (|4.13p has the formal solution 



^{t) =exp(-ii7oot)*(0) 

^[cos{fiBt) + 7^ sm{fiBt)]^{0). 



(4.14) 



Using Eqs. ([37T|) . ()4.12p and (|4.14p we then immedi- 
ately arrive to the following expression for the transition 
probability, P{t) — = sin^ {^Bt), which explains 

the behavior of the function P^™'''^-' at strong magnetic 
fields. Note that the analogous result was also obtained 
in Ref. 

Finally, it is worth of noticing that in contrast to the 
Dirac case, the behavior of the transition probability in 
the Majorana case is qualitatively similar for different 
matter densities and different magnetic fields [see Fig. 21 
panels (a), (c), and (e)]. 

The problem of Majorana neutrinos' spin-flavor oscil- 
lations was studied in Refs. [1, [9 [with help of numeri- 
cal codes. For example, in Ref. ^ the evolution equa- 
tion for three neutrino flavors propagating inside a pre- 
supernova star with zero metallicity, e.g., corresponding 
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X, 10^ km \iB, 10 ^'eV 



(e) (f) 

FIG. 4: (a), (c) and (e) The transition probability versus the distance passed by a neutrino beam. The neutrino parameters 
have the same values as in Figs. [l]|3l Eu = 10 MeV, Sm? = 8 x 10~^ eV'^ and 9 = 0.6. The magnetic energy is equal to 
= 10~^^ eV. The black line is the function P(x), the blue and green lines are the envelope functions Pu,d{x) and the red 
line is the averaged transition probability Pq{x). (b), (d) and (f) The dependence of the maximal values of the functions P{x) 
and Poix), blue and red lines respectively, on the magnetic energy for the given density. Panels (a) and (b) correspond to the 
matter density p — lOg/cc, (c) and (d) - to p = 50g/cc and (e) and (f) - to p = lOOg/cc. 



to W02Z model [2l[, was solved for the realistic mat- 
ter and magnetic field profiles. Although our analytical 
transition probability formula (|4.10p is valid only for the 
constant matter density and magnetic field strength, it 



is interesting to compare our results with the numeri- 
cal simulations of Ref. @ . In those calculations the au- 
thors used magnetic fields B ~ 10"'^'' G and magnetic 
moments ^ 10^^^ that give us the magnetic energy 
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liB ~ 10~ eV. This value is the maximal magnetic 
energy used in our work. 

It was found in Ref. Q that spin-flavor conversion 
is practically adiabatic for low-energy neutrinos corre- 
sponding to El, ^ 5 MeV inside the region where « 0.5 
and the averaged transition probability for the channel 
i>e is close 0.5. This big transition probability is 
due to the RSF-H and RSF-L resonances at the distance 
« 0.01 J?©. Even so that we study two neutrino oscil- 
lations scheme, we obtained the analogous behavior of 
pimax) ^ggg p.g^ ^ However, in our case this big tran- 
sition probability is due to the presence of the strong 
magnetic field (see Refs. [13, HI])- We cannot compare 
our transition probability formula (I4.10p with the results 
of Ref. [H for higher energies, E,, > 25 MeV, since spin- 
flavor oscillations become strongly nonadiabatic for these 
kinds of energies and one has to take into account the co- 
ordinate dependence of the matter density which should 
decrease with radius as [3^ . 

V. SUMMARY 

We have studied in this work the propagation of mas- 
sive flavor-mixed Dirac and Majorana particles in mat- 
ter under the influence of an external magnetic fleld in 
an approach based on relativistic quantum mechanics. 
The magnetic moment matrix is assumed to be both in 
the Dirac and Majorana case nondiagonal in the mass 
eigenstate basis. Starting with Lorentz invariant wave 
equations, the interactions of neutrinos with matter and 
magnetic fields included, we derived the effective Hamil- 
tonians and solve the evolution equations for the elec- 
trically neutral matter with n„ = Up. We then found 
the probabilities for the spin-flavor oscillation i/^ 
{v^ — > {v]^Y for Majorana neutrinos) and examined their 
dependence on the matter density and the strength of the 
magnetic field. 

We also studied in our approach the evolution of mixed 
massive Majorana neutrinos in vacuum, which is an ex- 
actly solvable problem. We derived the wave functions 
of the neutrinos and the transition probabilities for the 
process — > Va- We included the terms quadratic in 
neutrino masses and discussed the origin of these correc- 
tions of the standard formula for the probability. 

These corrections to the Pontecorvo formula, see 
Eq. p.l2p , are rapidly oscillating functions suppressed by 
the ratio ma/|k|, which is small for the relativistic neutri- 
nos. It was revealed in our previous works [Til, [l^. [Tsl [T^ 
that the analogous terms originate from the accurate ac- 
count of the Lorentz invariance in the study of oscilla- 
tions of Dirac neutrinos. It is also known that in de- 
scribing of evolution of Majorana neutrinos in vacuum 
the transitions of the types — > and {^aY 
are possible [33j . with the transition probability of the 
latter case being suppressed by the factor ma/|k|. The 
transitions — > {v]^Y can be interpreted as neutrino- 
to-antineutrino oscillations and manifest in the possible 



neutrinoless double beta decay [3J|. In our work on the 
basis of the relativistic quantum mechanics we derived 
the transition probability for neutrino-to-antineutrino os- 
cillations which is consistent with the results of Ref. [l^l ■ 

Our results can be used for analyzing the behavior of 
neutrinos in supernovae, in particular in the so-called 
zero metallicity presupernovae, where the condition n„ = 
Up is fulfilled [21] outside the Si-f O layer. It is known that 
these stars can possess very strong magnetic fields, up to 
10^^ G, or even stronger [35]. 

Wc found that large rates of transitions that change 
neutrino flavors and chiralities are possible both in the 
Dirac and Majorana cases with a typical oscillation 
lengths of the order of lO'^ — 10** km for matter densi- 
ties of the order of 10 — 100 g/cc, for typical supernova 
neutrino energies of 10 MeV and for (xB of the order of 
10^^"^ — 10~^^eV. The parameters of a shock wave, the 
density and width, can vary in a rather wide range, but 
typically they fit with the parameter values we used in 
our analysis [l^l- In particular, the width of the shock 
wave can be of about the same size as the oscillation 
length of the spin-flavor oscillations. It would thus be 
possible that the transition probability achieves its max- 
imum value (can be as high as 0.8 in some cases) when 
neutrinos pass through the shock wave on their way to- 
wards the outer layers of the star. If the flight distance 
were many oscillation lengths, the probability would be 
averaged to a smaller value. The results of our work 
can have the implication to the r-process nucleosynthesis 
since it is sensitive to the amount of neutrinos of various 
species emitted in a supernova [s^]- 

In the Dirac case the maximum value of the transition 
probability is achieved for a speciflc value of /^-B, while 
in the Majorana case the transition probability can be 
large for a wide range of /i_B values. Supposing that the 
neutrino magnetic moment is 3 x 10^ ^^/iB, allowed by 
the present astrophysical and cosmological data jsl] , the 
largest transition in the Dirac case would occur when 
the strength of the magnetic fleld is in the range 10® — 
10^ G. Although a magnetic fleld on the surface of a 
neutron star is typically stronger, 10^^ — 10"^^ G, at outer 
parts of the envelope it may have a suitable value for 
a large transition to take place. The difference between 
the dynamics of Dirac and Majorana particles under the 
influence of the same external fields can be a smoking 
gun to reveal the nature of neutrinos. 

Note that in the majority of the previous works de- 
voted to the spin-flavor oscillations the analytical tran- 
sition probability formulae were obtained only for Ma- 
jorana neutrinos. The case of Dirac particles was stud- 
ied only in connection with small off-diagonal magnetic 
moments. Our work studies the opposite situation of 
big off-diagonal elements of the magnetic moments ma- 
trix. The effect of the appearance of the big transi- 
tion probability at moderate magnetic field strength has 
never been described earlier. As for Majorana neutrinos 
the enhancement of the transition probability in strong 
magnetic fields is in agreement with the previous stud- 
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ies p^.[23|. Note that spin- flavor oscillations of Majorana 
neutrinos can be also resonantly enhanced in the moder- 
ate magnetic field, see Ref. Q. However this effect does 
not happen in isoscalar matter with ?ie = rtp = n„. The 
condition for the resonant spin- flavor precession reads , 
Gpine — «n) = ±(5m^ cos 20/4fc, which is not fulfilled in 
isoscalar matter. 

The obtained analytical formulae for the transition 
probabilities, Eqs. (|2:2T|) and (|4l0)) . are valid in the 
important case when Kg = 0.5, which is realized in 
the zero metallicity presupernova stars. These kinds of 
stars were quite common for early stages of the galaxies 
formation. We obtained that spin-flavor oscillations of 
Dirac neutrinos can be strongly influenced by the mod- 
erate magnetic field of a neutron star and the big tran- 
sition probability can exist not only for Majorana neu- 
trinos 0] . Thus spin- flavor oscillations of Dirac neu- 
trinos can significantly change the relic supernova neu- 
trino background 39]. Although now no signal of relic 
neutrinos was observed [40], there are still some efforts 



to calculate the flux of relic De for the KamLAND and 
Super-Kamiokande detectors [41|. The reliable simula- 
tion of the supernova explosion does not exist yet. For 
example, the propagation of a shock wave can signifi- 
cantly change magnitude and shape of the magnetic field 
which are very important in our calculations. Neverthe- 
less it is believed that the future galactic supernova neu- 
trino burst might give some clues to the physics of relic 
supernova neutrinos [1| . 
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